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Abstract — The maximum traffic arrival rate at tlie network for 
a given delay guarantee {delay constrained throughput) has been 
well studied for wired channels. However, few results are available 
for wireless channels, especially when multiple antennas are 
employed at the transmitter and receiver. In this work, we analyze 
the network delay constrained throughput of a multiple input 
multiple output (MIMO) wireless channel with time-varying 
spatial correlation. The MIMO channel is modeled via its virtual 
representation, where the individual spatial paths between the 
antenna pairs are Gilbert-Elliot channels. The whole system is 
then described by a A'-State Markov chain, where K depends 
upon the degree of freedom (DOF) of the channel. We prove 
that the DOF based modeling is indeed accurate. Furthermore, 
we study the impact of the delay requirements at the network 
layer, violation probability and the number of antennas on the 
throughput under different fading speeds and signal strength. 

Index Terms — Delay constrained throughput, correlated 
multiple-input-multiple-output (MIMO), Markov modeling, 
stochastic network calculus, moment generating function (MGF). 



I. Introduction 

Multiple input multiple output (MIMO) |]T| ||2) commu- 
nication has generated a lot of interest in recent years for 
providing high data rates and reliability without consuming 
extra power or bandwidth. It significantly improves the ca- 
pacity as compared to traditional single input single output 
(SISO) communication. The importance of MIMO technology 
can be well understood from the fact that IEEE standardization 
committee is actively pursuing IEEE 802.1 In and 802. 16e 
which promise high throughput using multiple antennas in 
wireless networks. 

Significant amount of work has been done on MIMO to 
predict the rate regions under different channel conditions, 
antenna settings, etc., under the assumption of infinite delay 
and queue / buffer sizes Q . Recently, however, there has been 
a rapid growth of delay sensitive applications like VoIP and IP 
video that are run on packet switched wireless networks. Albeit 
the decoding delay at the physical layer has been reduced by 
the invention of graph-like codes to the order of thousands of 
bits, the dela}[^at the higher layers has a significant impact on 
the overall performance for such systems. In addition, MIMO 
links have been considered for delay critical applications such 
as mobile ad hoc networks (MANETs) f3|, for they overcome 
the problems of low system throughput. It is therefore essential 

'Throughout the paper, the term delay refers to queuing delay and the 
throughput as the delay constrained throughput. 



to predict the rate region of a MIMO system for a given 
queuing delay guarantee. Indeed, the task of predicting the 
rate regions for a given delay guarantee is seen as one of the 
grand challenges in wireless networks |4|. 

In this paper, we formulate a method to find the delay 
constrained throughput of a spatial multiplexed correlated 
MIMO wireless channel on the flow level. We define the delay 
constrained throughput as the maximum packet arrival rate 
at the network for which the delay requirement is met. The 
service model is based on the Gilbert- Elliott (GE) channel |j5J, 
Q where for a given setting of the propagation environment, 
each channel path is considered in the virtual domain ||7], |[8| 
to be either in a good or a bad state. The good represents the 
case when there is a strong wireless link between the antenna 
pairs, while the bad corresponds to a scenario when the link is 
blocked by some large obstacle(s). The permutations of the Tx- 
Rx channel paths comprise the sub-states of the Markov chain 
describing the whole system. The sub-states are aggregated 
based on the degrees of freedom (DOF) available of the MIMO 
channel |j9| which results in a if-State Markov chain, where 
K depends upon the DOF. The DOF based modeling results 
in the reduction of the state space from exponential in the 
number of antennas to linear in the DOF. In addition to 
numerical simulations as in [j9j|, the method is justified here 
by a mathematical proof. 

Multi-antenna systems typically operate in one of the fol- 
lowing transmission modes. Spatial multiplexing is used to 
increase throughput by transmitting independent data flows 
over different antennas. Spatial diversity, on the other hand, 
increases reliability by transmitting redundant flows of infor- 
mation in parallel. There has been some work reported for 



MIMO diversity systems |[T0|, |[nj. In |[T0|, the concept of 
effective capacity has been used to carry out QoS analysis of 
independent and identically distributed (i.i.d) MIMO systems 
using finite state Markov chains to describe the received signal 
to noise ratio (SNR). However, it is well known ||7|, |[8| 
that the i.i.d fading model is highly idealistic — a realistic 
MIMO channel model is one in which the elements of the 
channel matrix are correlated. Furthermore, exploiting the 
DOF provided by the MIMO for spatial multiplexing is far 
more important for high-speed wireless systems than spatial 
diversity considered in [10|, [11]. 

A significant amount of work has been done on link-layer 



channel modeling of wireless channels p2) . In 1 13 1, diffusion 
approximation is used to carry out the throughput analysis of 



multihop memoryless wireless networks under SISO commu- 
nication. The concept of effective bandwidth is used in (14] 
to estimate the capacity of the wireless channels with memory 
that varies randomly with time and space. However none of 
these works covered the throughput analysis of more realistic 
spatially multiplexed correlated MIMO systems. 

The major contribution of this work is as follows. We 
present a methodology to calculate the network delay con- 
strained throughput of a correlated MIMO system that uses 
spatial multiplexing at the transmitter. In the analysis, we 
incorporate statistical independence of arriving traffic and 
the service provided by the channel. Secondly, we give a 
mathematical proof that the DOF based aggregation used in 
^ is indeed accurate. Furthermore we show the impact of 
increasing the delay guarantee and the number of antennas 
on the throughput under various conditions, such as signal 
strength and fading speed. 

The rest of the paper is organized as follows. We first 
present the correlated MIMO system model in Sec. |ll] A 
Markov modeling of the correlated MIMO system based on 



the DOF is carried out in Sec. Ill for calculating the delay 



constrained throughput. An application scenario based on 




IEEE 802. lln is presented in Sec. IV for which we calculate 
the numerical results. 

II. The System 

We consider a single user MIMO system with N transmit 
and M receive antennas and spatial multiplexing at the trans- 
mitter Throughout the paper, perfect channel knowledge is 
assumed to be available at the receiver while the transmitter 
has no information about the propagation environment. We use 
[•]^ to denote the matrix conjugate transpose, I„ denotes the 
identity matrix of order n, while Var( ) denotes the variance 
and CAf{n,tj}'^) denotes the circularly symmetric complex 
Gaussian distribution with mean /i and variance The Mx 1 
received signal vector y of a MIMO system reads 



(1) 



where x is a x 1 transmitted signal vector, n ~ CAf{0, Im) 
and p is the average total SNR at one receiving antenna. H 
is an M x N complex channel matrix with entries /i„i,n, 
where a particular hm,n represents the channel path from 
transmit antenna ji to receive antenna m. The channel gains 
are normalized so that E[|/im,np] = l,Vm,n. 

In most of the Uterature on MIMO channels (see for 
example (T] and fT\), the channel matrix H is assumed to 
have i.i.d complex Gaussian entries. The i.i.d MIMO channel 
model represents the ideal case where the environment is rich 
scattering and antenna spacing sufficiently large. It is known, 
however, that that this is seldom the case in practice |7|, |8 |. In 
this work, we use a more realistic model where the elements of 
the channel matrix are correlated and the correlation evolves 
in time according to a discrete time Markov chain. 

The virtual model proposed in ||7|, |j8| is a convenient 
method to link the physical propagation environment to the 



correlation structure of the MIMO channel when the corre- 
lation varies slowly compared to the Rayleigh fading experi- 
enced by the individual channel paths. The true channel under 
the virtual model reads 



H = AatHA 



t 



(2) 



where H represents the virtual channel related to H, while Ajy 
and Am are unitary discrete Fourier matrices (cf. |[7), |[8J). Of 
the many useful properties of H, the key observation is that the 
elements hm.n of H are independent random variables when 
the number of antennas is sufficiently large. Furthermore, 
the second-order statistics of /im,n are directly related to the 
scattering function in virtual domain and approximately given 
by the variance matrix V with entries Wm,n = Var(/i,„ „). 
Since the original channel matrix H is normalized to have 
unit power entries, it follows from (|2| that the variance matrix 
must satisfy 



NM 



(3) 



It is to be highlighted that the virtual channel H with a 
variance matrix V generates a correlation structure on the true 
channel H as given in fTS, Eq. (9)]. This in turn maps to a 
specific physical scattering environment. For simplicity, we let 
the scattering environment and hence the spatial correlation 
matrix, change between code words and consider the fading 
experienced by each of the code words to be ergodic. This 
models a scenario where the positions of large objects (e.g. 
buildings) between the user terminal and the base station 
change relatively slowly, whereas the scattering environment 
around the user (induced e.g. by cars) experiences rapid 
changes. If we let the MIMO channel H be Rayleigh fading, 
the virtual channel matrix H has (approximately) independent 
complex Gaussian entries hm,n ^ CJ\f{0,Vm,n)- Since we 
assume no CSI at the transmitter, the inputs are chosen to 
be independent standard complex Gaussian and the ergodic 
capacity conditioned on V reads (S) 



C(V) = E log 



det ( Im + ^HH1 



(4) 



where the expectation is with respect to all Gaussian channel 
matrices H that have variance matrix V. 

III. Delay Constrained Throughput Analysis 

The capacity ^ depends on the scattering environment via 
matrix V. Instead of considering one specific V for all trans- 
missions or picking it e.g. randomly, we let the propagation 
environment evolve between the code words according to a 
discrete time Markov chain with state space S. Note that 
although we focus in the following on the simplified setting 
where the number of antennas at both ends is N, generalization 
to arbitrary M x N MIMO channels is also possible. We shall 
omit these cases for the ease of exposition. 



A. Markov Modeling 

The Markov chain describing the scattering environment is 
defined through the individual spatial paths by modeling them 
as identical and independent Gilbert-Elliott (GE) channels in 
the virtual domain. We use the notation S{hm.n) to denote the 
state of the link (to, n) and consider for simplicity only two 
states, good ('g') or bad ('b'), so that 



(5) 



S{hrn,n) G {b,g}, Vto, n = 1, 

The sM^J-^stafe^s e 5 of the Markov chain describing the time 
evolution of the variance matrix V are then written as 



s= ••• SihN^N)] e {b,g} 



(6) 



For the GE channels, the good state corresponds to a strong 
wireless link between the given antenna pairs (in the virtual 
domain), while a bad channel represents the case when the link 
is blocked by some large obstacle(s). For simplicity, we follow 
the framework of p5) and let the bad state of the link (to, n) 
be defined as Vm,n = 0, while the good state corresponds to 
Vm.n — Vs> I- Note that the variable Vg depends on the sub- 
state s e 5 of the Markov chain but not on the path indices, 
and is selected for all s G 5 so that ^ is satisfied. Albeit this 
is the simplest possible scenario and selected for expository 
reasons, extensions to non-binary V are possible. 

The probability for a particular link (to, n) to go from the 
good state to the bad state is given by Pgh and from the bad 
state to the good state by pbg- For a certain path, the steady 
state probability for the bad channel state reads thus 

n=^^. (7) 
Pgb + Pbg 

It is to be noted that the Markovian model of the time varying 
correlated MIMO channel is completely specified by the two 
independent parameters k and pbg which describe the fading 
speed, i.e., the speed of change of the scattering environment. 
These parameters in turn can be adjusted for example to meet 
the results of a measuring campaign. 

B. State Aggregation Based on DOF 

Since all individual links (to, n) are independent two-state 
GE channels in the virtual domain, there is by definition a 
total of \S\ ~ 2^ sub-states in the Markov chain describing 
the scattering environment. We call this method : 1. The 
problem of this approach is, however, that the state space 
grows exponentially with the number of spatial paths iV^, 
inducing a high computational complexity in simulations. To 
alleviate the problem, we consider next an approximation for 
method : 1 that reduces the state space significantly. 

In what we call method : 2, the set of sub-states S is divided 
into K disjoint subsets {Si \ i — 1, . . . ,K}. Each sub-state 
s G Si corresponds to a variance matrix V that provides r = 
i — 1 spatial DOF. This leads to a new Markov chain with 
states i = 1, . . . , K — N + 1, providing a significant reduction 

-We use the term sub-state instead of state in order to separate tlie state 
space of two different Markov cliains. Details are given below. 




Fig. 1. States of the Markov Chain for Af = 3 



in the state space compared to method : 1. We remark that 
method : 2 was used to our knowledge first time in this context 
in |[9), where it was found to give accurate results via extensive 
numerical simulations. In this paper a theoretical justification 
for method : 2 is provided by the following lemma. 

Lemma 3.1: Let V be an arbitrary N x N binary variance 
matrix for which the degrees of freedom are given by 



r(V) 



1 = rank(H), 



2,3, 



,N+1, (8) 



with probability one when the power constraint (|3]l is pre- 
served. Then, for sufficiently large r(V), the ergodic capacity 
at high SNR is given by 



C(V) = r(V) 



log2 



pN 



c'(V) 



(9) 



where c'(V) e 0(1) is a constant that does not depend on 
SNR, and r(V) corresponds to the degrees of freedom, or 
multiplexing gain, of the virtual channel H whose second 
order statistics are given by the variance matrix V. 

Proof: The proof is deferred to the Appendix. ■ 

Remark 1: Lemma 13.11 is an extension of the well known 
high SNR capacity result [ 1 1 where Wm,n = 1 Vto, n, to the 
case of arbitrary Vm,n G {0, Vs]^ s e 5. As a consequence, for 
the state i in method : 2, the ergodic capacity of the system 
grows cx (i — l)log2(p), while the differences between the 
sub-states are vertical shifts that become negligible as p cxi. 

Fig. [T] depicts the ergodic capacity of a 3 x 3 MIMO system 
and shows that the conclusions of the Lemma ITT] are accurate 
also for small r. As noted in Remark[T] the difference between 
the sub-states manifests as a constant vertical shift at high 
SNR. The thick lines in Fig. [T] represent in fact bundles of 
sub-states. 

Definition 3.2: Let Si be the set of sub-states and \Si\ its 
cardinality for the state i = 1, . . . , K in method : 2. We define 
the transmission rate Ri in state i as 



Hi — 



1 

]S~ 



Vi = 1, 



,K, 



(10) 



seSi 



where Cg is the average capacity of the sub-state s ^ Si given 
by (|4]). The set of all rates {Ri | i = 1, . . . , K} are collected 



in the diagonal rate matrix R(0) = diag(e 



)■ 



Definition 3.3: Let Eij be the event that corresponds to the 
transition of the Markov chain from state i to state j and let 
Qij be the corresponding transition probability 

J - 51 • J I ^) ' Vz,j = l,2...if , (11) 

seSi 

where the normalized steady state probability tTs of being in 
sub-state s is given as 



(12) 



and 



N 

while for all ijj — 1^2 ... K 



e{K,i-K} , (13) 



N 



Pr \ ses,)= J2 n P-y ^ ^i^) 

s'^Sj m,n— 1 

where k denotes the element of the vector and pr^y E {pgb , 1 — 
Pgb,Pbg, 1 — Phg}- Note that the entries qij in ( [TT| comprise 
the elements of the transition probability matrix Q. The steady 
state vector tt = [tti tt2 ■ ■ ■ ttk] can then be obtained by 
solving TT — ttQ, and we define the first-order capacity 
as 



K 



(15) 



1=1 



Example (2x2 MIMO system): Let us consider a MIMO 
system that has 2x2 antenna setup and uses spatial multiplex- 
ing at the transmitter. The correlation structure of the channel 
matrix H is induced by the virtual domain variance matrix V 



and can be obtained from 1 15 Eq. (9)], as explained earlier 
Here, however, we are not interested in the specific forms 
of spatial correlation, but rather describe how to model this 
simple MIMO system by a three state Markov chain for which 
the states 1,2 and 3 correspond to a transmission scenario 
where the channel has DOF of 0, 1 and 2, respectively. 

Let us write the sub-states in method : 1 explicitly in the 
form s = S'(/ii,2) S{h2,i) S'(/i2,2)]- The method : 2 

groups the sub-states s e 5 of the original Markov chain as: 

. State 1: [bbbb]. 

. State 2: [bbgg], [bgbg], [gbbb], [bgbb], [bbgb], [bbbg], 
[gbgb], [ggbb]. 

. State 3: [bggg], [gbbg], [bggb], [gbgg], [ggbg], [gggb], 
[gggg]- 

As we saw in Lemma |3.1| and Fig. [T] the rate at the physical 
layer is approximately equal for all sub-states that fall within 
the same DOF when V is a binary matrix. The transition 
matrix Q of this three state process can be calculated using 



([TT]) - ( [T4] i. The steady state vector follows accordingly and 
is given by 

(1 - kY + 4k(1 - nf + 2k2(i - i^f] e rIxs . 

The corresponding rate matrix reads 

R(0)=diag(l,e''«%e''«^), 

where we assume that there is no workload processed in 
state 1. 

C. Delay Constrained Throughput 

We are now in a position to present the main result of 
this paper The moment generating functions (MGF) of a 
stationary process X{t) is given by Mx{0,t) — E [e'''^'*)], 
where E [•] is the expectation. In the sequel, we denote 
Mx{0,t) = Mx{-d,t) for pai-ameter 6* > 0, t > 0. The 
discrete time arrivals and service of the channel are assumed 
to be statistically independent and stationary random processes 
and are given by the cumulative processes A{{),t) and S{s,t) 
respectively. For all real 9, the arrivals and the service random 
process posses the MGFs M^(0,i) and Ms{0,t) respectively. 

Lemma 3.4: The MGF of the random process S{t) de- 
scribed by a homogeneous Markov chain, with transition 
matrix Q and stationary state distribution vector tt, is given 
for 6* > and t > by 

Ms{d, t) = 7r(R(-0)Q)*-iR(-0)l, (16) 

where 1 is column vector of ones while 

R(0) =diag(e''«S...,e^-«^). (17) 

Proof: We refer the interested reader to p6) for the proof. 



We denote the traffic arrival rate by A while d\ represents a 
bound with violation probability e on the delay. We can only 
provide a delay guarantee d'^ if d\<d'~^ . 

Proposition 3.5: The delay constrained throughput Ad of a 
correlated MIMO system under a delay guarantee d'^ is given 
as 

Ad = max {\\dl< d^}, (18) 
where d\ assuming FIFO scheduling is given as 



= inf 

9>0 



inf 




UA{9,s-T)Us{B,s)-\TLe\ <0 



(19) 



and \\l\s{9,s) is given by ( fTSj l. 

Proof: We refer the interested reader to p6| for the proof 
of^. ■ 
Remark 2: The delay bound ([T9]l is calculated using 
stochastic NetCal approach based on MGF fTTl. MGF for 
a variety of arrival models is available in the literature | [T6[ 
but finding the MGF of the service process \\l\s{9^t) is a 
challenging task. We address this challenge by making use of 



Lemma 3.4 to calculate Ms{9,t). Having obtained Ms{0,t), 
a stochastic bound on the delay d\ in ( [T9] l can be obtained 
using Chemoff's bound, Boole's inequality and applying the 
technique in ||T7|. 

We next use the result in this section to calculate the delay 
constrained throughput of a x MIMO channel which is 
parameterized according to IEEE 802.1 In. 

IV. Application and Numerical Results 

We present numerical results for the delay constrained 
throughput of a MIMO system which is parameterized ac- 
cording to 802.1 In using 40 MHz channel |18|. We choose 
a base time unit of 31fis for the discrete time model which 
is approximately the time needed to transmit a data block of 
2312 bytes over a 600 Mbps channel. The rate matrix R(0) 
is obtained through an extensive simulation as described in 



Section III-B We omit showing confidence intervals as these 
turn out to be very small. We fix the SNR to 25dB for all 
N and set the rate of the highest state rx corresponding 
to this fixed SNR. The rate of the rest of the states follows 
accordingly corresponding to the different channel matrixes. 
The infinite sum in the delay bound formula ( [T9] l is computed 
for the first 4000 units of time. For the GE model, we set 
the transition probability from bad to good pbg = 0.1 and the 
probability from good to bad — 0.01 similar to [19|. 

The framework discussed earlier enables the derivation of 
delay constrained throughput of a various number of traffic 
sources whenever the MGF exists. A collection of such MGFs 
for sources of different statistical properties can be found in 
In this work we use a periodic source for the discrete 
time model that generates arrival traffic. Such a traffic source 
with period t produces a units of workload (data blocks) at 
times {Ut + tit, n = 0, 1, . . .} where U is the starting time 
which is uniformly distributed in the interval [0, 1]. The MGF 
of such a source is known (see e.g. pO)) as 



MAiO,t)=e' 



1 



r 



(20) 



for t > and 9 > 0. We parameterize the source such that 
its period r is 10 time units and set the number of generated 
data blocks a to achieve a given arrival rate A. 

Fig. |2] depicts throughput as a function of delay guaran- 
tee for different delay bound violation probabilities e 
using method : 1 without DOE aggregation (WODOF) and 
method : 2 with DOF aggregation (WDOF) for = 2 
antennas. It shows that the method : 2 of aggregating the 
states based on DOF is reasonably accurate. For the results 
to follow, we use method : 2 of DOF based aggregation only. 
Furthermore, one can observe the decrease in throughput as 
the violation probability gets tighter. 

It is to be noted that Shannon limit is the upper bound to the 
delay constrained throughput as can be seen in Fig. |3]for N = 
2, where we fix e = 10^^. Also note that the improvement 
in throughput is less as we go to the higher values of delay 
guarantee. 




100 150 200 

delay guarantee [time slots] 



Fig. 2. Delay constrained throughput for different violation probabilities e 
using method : 1 (WODOF) and method : 2(WDOF) 
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Fig. 3. Delay constrained throughput for different SNR's 



Tab. |l] verifies the known notion that higher antenna config- 
uration provides a higher capacity. It is to be noted that the 
ratio of the first-order capacities of = 3 to = 2 stays 
unchanged over different fading speeds. 

We next study the effect of fading speed on the throughput 
for different number of antennas (A^ — 2, 3) in Fig. |4] on 
loglog scale where the delay guarantee is fixed to 30 time 
slots and e = 10"'^. We fix the steady-state error rate k for 
each spatial path given by (|7]i to 0.1. We control the fading 
speed, under fixed k, by the mean time for the Markov chain 
of the path to change states twice. To this end we variate pbg 
such that a small value corresponds to slow fading and vice 
versa. The notable thing is that for slow fading the difference in 
throughput between N = 2 and A^ = 3 is large and for higher 
fading speeds, which corresponds to a channel nearly without 
memory, the profit from extra spatial paths diminishes. As the 



TABLE I 

First-order capacity for different # of antennas 



N 


[bps/Hz] 


2 


14.1 


3 


20.82 




10"' 
fading speed p 



Fig. 4. Throughput as a function of fading speed 
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Fig. 5. Throughput as a function of different # of antennas A'^ 

first order capacity given by ( fTS) and shown in Tab. |l] stays 
constant over the range of fading speeds, which is implied by 
the fixed k, the difference in throughput might be due to the 
higher order statistics of capacity. Furthermore, the impact of 
statistical multiplexing between the spatial paths, which is 
higher for = 3, is apparent in the slow fading regime where 
the Markov chain changes states less frequently. 

It is well known fT\ that an increase in the number of 
antennas leads in the best case to a linear increase in channel 
capacity. This increase capacity manifests itself in the increase 
of the delay constrained throughput for a fixed scenario 
conditions. This can be seen in Fig. [s] where e — 10^^. 
It saturates after a given delay guarantee owing to the well 
understood concept that traffic arrival rate should never be 
greater than the capacity to ensure system stability. 

V. Conclusions and Future Work 

In this paper we formulate a method to find the delay 
constrained throughput of a correlated, spatially multiplexed 
MIMO system in which the MIMO wireless channel is mod- 
eled as a Markov chain. The states of the Markov chain 

^Statistical multiplexing, which is a flow level phenomenon, should not to 
be confused with spatial multiplexing of the MIMO system. 



are decided on the degree of freedom for which we provide 
a rigorous mathematical proof. We parameterize the model 
to the application of IEEE 802. 11 n. We present numerical 
results where we quantify the impact of increasing the delay 
guarantee and the number of antennas on the delay constrained 
throughput for varying signal strength and fading speed. 

This work finds it application in performance evaluation 
of wireless networks where we are interested in the maxi- 
mum attainable throughput for a given delay guarantee and 
resources of the network. While in this work, we only used 
periodic source to illustrate the results, the same methodology 
can be applied to any traffic source with known MGF The 
formulation can also be extended to multihop and multi-user 
scenarios and it is left as our future work. 
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Appendix 



Here we sketch the proof of Lemma 3.1 introduced in 



Section |III-A| The result is obtained in the large system limit 
where the (effective) number of antennas on both sides grow 
without bound. We consider only N x N Gaussian MIMO 
channels that are described in the virtual domain by a binary 
variance matrix V with r non-zero rows and columns. The 
derivation for general M x N MIMO channel as well as the 
case of unequal number of non-zero rows and columns follows 
by the same procedures but is omitted in the present paper due 
to space constraints. We also remark that our result is akin to 
the one proposed in (15], but the proof is slightly different 
due to the special structure of the matrix V considered in this 
paper. 

Lemma 5.1: Let r = rank(H) with probability one, and let 
V be such an N x N binary variance matrix that by deleting 
the all-zero rows and columns from it we get a modified r xr 
variance matrix V. Denote the number of non-zero elements 
in the row m = 1, 2, . . . , r of V by D„i, let D = X]m=i 
be the total number of non-zero elements in V (and V) and 
write V — {Dm ■ m — l,...,r}. Then, for fixed V, the 
ergodic capacity reads 



log^ipN) - 
1 

^ rln(2) 



1 



7 



ln(2) 

r 



(21) 



ciV) 



where r — ?► oo, 7 = lim„^co {Hn - ln(n)) = 0.57721566 . . . 
is the Euler-Mascheroni constant and iJ„ = J2p=i 
nth harmonic number. 

Proof: Let H e cxr ^ modified Gaussian virtual 
MIMO channel that is obtained from H by removing the all- 
zero rows and columns. The non-zero elements > : 
m,n = 1, . . . , r} of H are independent by assumption |8 1 and 
their variances are given by the corresponding elements of V. 



The ergodic capacity per DOF at high SNR, given V, can then 
be written as ||T5| Eq.(78)] 

^ = J Elog, [det (HHt)] + log, (^) + O ( A) . 

(22) 

Note that the last term in ( p2j i vanishes in the high-SNR limit 
p — ^ oo. The second term depends on the total number of 
transmit antennas N since the transmitter has no (statistical) 
channel information and uses therefore uniform power alloca- 
tion over all antennas. By |15j Lemma 5], 



lim -Elog2 [det (HH^)] = lim - j logs 



~D 



rln(2) ^ 

^ ' m=l 




where xio ^ chi-square random variable with 2Dm 
degrees of freedom. The expectation in ( |23] l has a closed form 
solution (see e.g., pT| Lemma 10.1]) 



so that plugging (|23 
high SNR 



EMx^Z5„J-^i3„.-i-7 , (24) 
'24\ to p2l) and simplifying yields at 



lim^ 

r— J-oo T 



lim < log2 



pN 

1j 



Hn(2) ^ ^^--^ f 

^ ' m=l ) 



1+7 

(25) 



completing the proof. ■ 
Remark 3: For large degrees of freedom r, and a square 
(modified) variance matrix V, the ergodic capacity (|4]i depends 
on V only through c(I?) in pT| ). The exact structure of V is 
thus not important, just the number of non-zero paths between 
the receive antennas and the transmit antennas. Furthermore, 
for fixed r and V, the maximum capacity difference between 
the sub-states can be explicitly calculated from c(I?). In the 
following, however, we do not consider the exact calculation of 
c(I?), but rather study its behavior in the high-SNR and large 
system limits. This will provide and expression for capacity 
that holds for all choices of T> as stated in Lemma 13.11 



Proof of Lemma 3.1 



Let us first consider the case when V is not sparse, that is, 
e 0{r) or 



Dm 
iim 

r— f oo T 



< Um < 1, 



m = 1, 



,r . (26) 



From Lemma |5.1| and the bound p2| 
1 



24(n 



-yp <i/„-(ln(n + l/2)+7) < 



24^2 



(27) 



we have 

lim c{V) = 

r— J-OO 

+ lim 



ln(2) 



1 



•ln(2) ^ 

^ ' m=l 



In 



1/2 



D 



O 



D2 



(28) 



where the last term in ( |28] l vanishes rapidly with {DmYm=i- 
By definition 

A„ - 1/2 



lim 
SO that 

lim c{V) = 



D 



lim 



Er J 
n=l ^ri 



7 
ln(2) 



lim - 

m— 1 



l0g2 



(23) Plugging d29b to Lemma 5.1 shows that as r 



C{V) 



1 



1 



ln(2) 



l0g2 



m— 1 



an 
oo, 

a. 



(29) 



(30) 

Thus, for all V that satisfy ^ and (|26]l, at high SNR and for 
large but fixed degrees of freedom r, the ergodic capacity is 
given by 



CiV) = r 



log2 



pN 



+ c'{V) 



(31) 



where c'{V) e 0(1). 

By similar calculus as above, it is not difficult to show that 
pT] ) holds also if the matrix V is sparse, that is. Dm € o(r) 
or 

lim — = 0, Vm = l,...,r , (32) 

r— >oo T 

where the I?m's are finite constants or grow sub-linearly with 
r (for example, = ^/r). Furthermore, Lemma 



5.1 



can be 

modified for non-square V, and the rest of the proof follows 



a similar calculus as shown here. The final result is that (31 



holds also for such variance matrices, as stated in Lemma |3TT| 
These cases lead, however, to more cumbersome notation and 
we have omitted them in the proof due to space constraints. 
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